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GENERALIZATION OF SOME DETERMIN ANTAL IDENTITIES FOR 
NON-SQUARE MATRICES BASED ON RADIO’S DEFINITION 

A. AMIRI \ M. FATHY \ M. BAYAT 1 2 


Abstract. In this paper, we focus on Radio’s definition for the determinant of non-square 
matrices. We develop some important properties of this determinant. We generalize several 
classical important determinant identities, including Dodgson’s condensation, Cauchy-Binet, 
and Trahan for non-square matrices. Also, we propose an efficient algorithm with 0((mn) 2 ) 
time complexity for computing Radio’s determinant based on Dodgson algorithms and dynamic 
programming technique. 
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1. Introduction 

The determinant is a fundamental concept in linear algebra and applied statistics, and has 
significant applications in various branches of mathematics and engineering. Originally deter- 
minants are defined only for square matrices. Therefore, if one can generalize the definition 
to non-square matrices, then, we will be able to extract some of the important properties and 
identities of non-square matrices that can be used in defining many other important concepts 
such as eigenvalues, eigenvectors and efficient matrix decompositions for these matrices. 

To date, many definitions have been proposed for the determinant of non-square matrices. 
Earlier works have been mainly focused on utilizing the determinant of square blocks to define the 
determinant of the non-square matrix. Joshi [6] defined the determinant of non-square matrices 
as the sum of the determinants of square blocks. They studied many useful properties of this 
determinant. Radic [9] proposed an efficient definition that has some of the major properties of 
the determinants of square matrices. Also, some other properties of Radio’s determinant and its 
geometrical interpretations, involving polygons in the plan R 2 and polyhedra in R 3 are given in 
[8]- [13] . Also, Amiri et al. [1] defined generalized eigenvalue and eigenvector based on Radio’s 
determinant. They applied it on video retrieval and video shot boundary detection applications 
[1]- [3] . Very recently, Y. Haruo and T. Yoshio [5] defined non-negative determinants for non- 
square matrices that have many valuable applications in multivariate analysis. They studied a 
broad spectrum of properties of this determinant. 

Successful as they are, existing works on the determinants of non-square matrices have not 
focused on generalizing many useful square determinantal identities, such as Dodgson, Cauchy- 
Binet and Trahan, for non-square matirces. Also, there is no efficient algorithm for computing 
these determinants. With this in mind, in this paper we focus on Radio’s definition [9] and 
we set out to develop some important properties of this determinant and to generalize some 
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of important determinant al identities for non-square matrices based on this definition. We will 
develop an efficient algorithm for computing this determinant based on generalized Dodgson 
identity and a dynamic programing technique that has polynomial time complexity. 

The remainder of this paper is organized as follows. Section 2 reviews Radio’s definition for 
determinant of non-square matrices. Section 3 studies some properties of this determinant. We 
extend the Dodgson identity for non-square matrices. We also design an efficient and straight- 
forward algorithm for calculating this determinant with polynomial order time complexity in 
Section 4. Section 5 presents the generalized Cauchy-Binet and Trahan formulas for non-square 
matrices. Section 6 summarizes the paper. 


2. Radio’s definition 


In the following, we first present Radio’s definition [9] for determinant of non-square matrices. 
Definition 2.1. Let A = (o n ,j) be an m x n matrix with rn < n. The determinant of A, is 
defined as: 


det(A) = J2 (— l) r+s det 
l<jl<—<jm<n 


aiji 


Clin j i 



a mjm 


(1) 


where ji, j 2 , • • • ,j m £ N, r = 1 + 2 + ■ ■ ■ + m and s = j\ + • • • + j m . If m > n, then we define 
det (A) = 0. 

Example 2.2. 


det 


Ol a 2 
bi b 2 


C13 

bs 


= (- l) (1+2)+(1+2) det 


ai a 2 
b± b 2 


+ (-l) (1+2)+(1+3) det 


a i 03 
bi 63 


+ 


+ (- l) (1+2)+(2+3) det 


02 03 

b 2 63 


= det 


ai a 2 
bi b 2 


det 


01 03 

h b 3 


+ det 


02 03 

b 2 63 


The determinant (1) is a skew-symmetric multilinear functional with respect to rows (only 
multiaffine with respect to columns) and, therefore, it has many well-known standard properties. 
The main properties are listed in the following theorem: 

Theorem 2.3. Suppose that A is an m x n matrix with m < n, then 

(1) If A = [ai,a 2 , • • • ,a n ], then det(T) = ai — a 2 + . . . + (-1 ) n ~ l a n . 

(2) If a row of a matrix A is multiplied by A, then the determinant of the new matrix is 
A det(yl). 

(3) If the matrix A has two identical rows, then det(vl) = 0. 

(4) If a row of matrix A is a linear combination of some other rows, then det(Al) = 0. 

(5) The interchange of any two rows in A will only change the sign of det(vl). 

(6) The Laplace expansion theorem is valid. 

Proof. See [1-3], [8], [10], [12], 

In addition, some of the other properties of Radio’s determinant and its several excellent 
geometrical interpretations involving polygons in the plan R 2 and polyhedra in R 3 are studied 
in [8]-[10-12], 


3. Other properties of Radio’s determinant 


The following lemmas and theorems exhibit some of the important properties of the Radio’s 
determinant. They are essential in studying of the theorems of the subsequent sections. 
Lemma 3.1. Let A be an m x n matrix, 1 < m < n, and m + n odd, then: 


det 

1 <i<m 
1 <j<n 


(a 




= 0, 
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where aij = 1 for all j, 1 < j < n. 

Proof. We proceed by induction on even integer n for all odd integer numbers m, 1 < m < n. If 
n = 2 then m = 1 , we have 

det ([ 1 , 1 ]) = 1-1 = 0 . 

Now, we assume that it is true for even n and all odd m, 1 < m < n. We will show that the 
identity holds for n + 2 (which is even number) and all odd m, 1 < m < n + 2. By expanding 
the determinant with respect to row m, we get 




n+2 


r -1 

det 

l<i<m 


= ^(-l) fc+m a m , fc 

det 

l<i<m 

(Oj j) i+7Tl 

j + k 

1 <j <n+2 



l<J<n+2 



Next, by expanding the right determinant with respect to the row m — 1, we obtain 


det 

1 <i<m 
1 <j <n+2 


j+k 


71+2 


= ) J"!) 


k'+m—l 


k= 1 
k'^k 


®m- i,k' det 

1 <i<m 
l<j<n+2 


(CLi.j ) ij£m — l,rr. 
3^k,k' 


therefore, we get 


det 

1 <i<m 
l<j<n+2 




n+2 n+2 

-tFi-ir 1 


k = 1 v=i 
fc+fc 


Q j m,kQ j m— l,k' 


det 

l<i<m 
l<J<n+2 L 


( CL'i j ) ij£m—l,rr 
j^k,k' 


On the other hand, all determinants on the right hand side of the above equality are of size m . 
by n, with the first row entries equal to one. Thus, by induction hypothesis, those determinants 
are zero and the proof is completed. Another case can be established in the same way. 

Lemma 3.2. 

det [Ai, ■ ■ • , A n ,0 m ] = det [+, ■ ■ • , A n ], 

and, 


det [A i , • • • , Aj—i, 0 )7l , Aj+i , • • • , A n ] — det [A \ , • • • , Aj —\ , Aj+i , • • • , A n ] , 
where m < n, A\ z = [ai^, • • • , a m ^] T for k G {1, ■ • • , n} — {j}, and 0 m is an m by 1 zero vector. 


Proof. The proof of the first formula is straightforward by Definition 1.1. We prove the second 
formula by induction on n. Set 


+ = 


B, 


, 


B: = 


°2 ,j 


‘mj 


(1 < j < n). 


Expanding the determinant with respect to the first row yields 


D — det [Ai, 5 Aj—i) 0 m , Aj- |_i, • • • , An\ 


det 


0 * 1,1 • • • CL\j — i 0 nij+i * * * 0*1,71 
B\ • • • Bj— i 0 m _i Bj+i • • • -B n 


— ^ 1,1 det [-B 2 , 5 -Sjf — 1 5 O 777, — 1 5 -®j+i ? * * * * -®n] * H - 

+(— l) J_ 1 + 1 aij-i det [Si, • • ■ ,Bj_2,0 m _i,Sj + i, • • • . B„\ - 
+(-l) J ’ + 1 + 1 aij + i det [Bi,-- 1 Bj— i, 0 m _i, Bj. |_2, , B n ] + • • • + 

+(-l) n+ 1 ai in det [Bi, • • • ,Bj-i,0 m -i,B j+ i,--- ,B n - 1 ] , 
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next, by induction hypothesis 


D — c i . ] det \B 2 , , Sj—i, , B n ] T 

+(-l)- 7_1+1 aij_i det [Bi, • • • , Bj- 2 , ~B j+1 , • • • , -B n \ + 

+(— l)- ,+1+1 aij+i det [B\, • • • , Bj_ i, —Bj + 2 , • • • , —B n ] + • • ■ + 
+(— l) n+1 ai,ndet [B \ , • • • , Bj 1 . — />’j . 1 . • • • , — -B n -i] 


or, equivalently 


.D = det 


Ol,l 


«ij-i — oij+i 
— -Bj+i 


— det [A\j • • • , j4j_i, , d. n ] . 


&gn 

--Bn 


Theorem 3.3. Suppose 1 < m < n, and m + n be an odd integer, A = ( djj ) 
an m by n matrix, and X be an arbitrary m by 1 column vector, then: 


det [A\ + X, ■ ■ ■ , A n + X ] — det [A\, • • • , A n ] . 


[^4i, - • - ,A n \ be 


Proof. We will prove the assertion by mathematical induction on n even and all m odd 1 < 
m < n. For n = 2, m = 1, we have 


det [dpi + x, a ,\.2 + x] = det [dpi, ap 2 ] = dpi - dp 2 . 


Now, assume the assertion is true for all n even and m odd with 1 < m < n. We prove it is also 
true for n + 2 even and all m odd with 1 < m < n + 2. By setting 


Aj = 


Bi 



a m,j 


(1 < j <n + 2) 


we get 




X2 

X\ 

X' 

, x' = 




Xm 


— det [A\ + X , • • • , A n + 2 + -T] — det 


cq.i + aq • • • dp n + 2 + X\ 
^Bx+X' ••• Bn + 2 + x' 

= (-1) 1+1 (dpi + xi) det [S 2 + X', ■ ■ • , B n+ 2 + X'] + ■ ■ • + 

+ (-l) 1+n+2 (op n+2 + aq) det \Bi + X',-- - , B n+1 + X'} , 


now, by inductive hypothesis 

D = (-l) 1+1 (dpi + xi) det [B 2 , • • • , B n+2 ] + 

+ • • • + (— l) 1+ " +2 (di in _|_2 + aq) det [B\, ■ ■ ■ , -B n+ i] = 

= {(-1) 1+1 (dpi) det [B 2 , ■■■ , B n+2 \ + • • • + (-l) 1+n+2 (op n+2 ) det [B u ■ ■ ■ , B n+1 ]} - 
— ci |(— 1) 1+1 det [B 2 ,--- , B n+2 \ + • • • + (— l) 1+n+2 det [B\, • • • , -Brj+i] } = 


= det 


opi 

B\ 


O-l, n+2 

B n 


X\ 


1 

Bi 


1 

B n +2 


considering Lemma 3.1, the last determinant is zero and the proof is completed. The second 
case can be treated similarly. 

Corollary 3.4. For m + n odd, 1 < m < n and for all k, 1 < k < n, we have 


det [A \ , • • • , A^_i, Afc, Afc+x , • • • , A n ] det \A k A \ , 


, Au — A 


k— 1 ? 


^4fc+i — A k 


1 A n Afc\ . 
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Proof . Applying Theorem 3.3 with X = — A*, and Lemma 3.2 the result is obtained. 

Theorem 3.5. [Cyclic] If 1 <m < n, and m + n is an odd integer, then for all i &{!,■■ ■ , n}, 
we have 

(_l)(i+i)m det[i4i) ... ,A i . 1 ] = det[A u — ,A n ], 

Proof. It is sufficient to prove 

(— l) m det[A n , Ai, ■ ■ ■ , Ai-i] = det[Ti, • • • ,A n \. 

Applying Theorem 3.3 with X = —A n and Lemma 3.2 we have 

(— l) m det[A„, Ai, • • ■ , A n _i] = (— l) m det[O m , A\ — A n , • • • , A n -i — A n \ = 

= (-ir(-l) m det[Ai - A n , ■ ■ ■ , A n -\ - A n ] = 

= det [A\ - A n , ■ ■ ■ , An-t - An, A n - A n ] = det[Ai, • • • ,A n ]. 

Theorem 3.6. [Senri-Cyclic] If 1 < m < n, and m + n is even, then for all i £ {1, • • • , n} 

( i)( n— ®) m det[j 4j, ,A n , —Ai, ■■ ■ , —Ai-i\ = det[A l5 • • • , A n ]. 

Proof. It is sufficient to prove that 

det [A n , —A\, ■■■ , -A„_i] = det [A\, ■■■ , A n ). 

Applying Theorem 3.3 with X = — A n and Lemma 3.2 we have 

det [An, Ai, , A n _i] — det[A n , A\,--- , A n —\, 0 m ] — 

= (— l) m det[— Ai, • • • , — A n _i,O m , A n ] = 

= (— l) m det[— Ai, • • • , -A n _i, -A n ,0 m \ = 

= (— l) m (— l) m det[Ai, ■ • • ,A n ] = det [A \ , • • • ,A n ]. 

The next theorem involves the product property of the determinant of non-square matrices. 
Theorem 3.7. Let 1 < m < n, and A be an m x m matrix, and B be an m x n matrix, then 

det (A.B) = det (A). det(-B). 

Proof. Let B = [B \ , • • • , B n ] , then 

det (A.B) = det ^A.[B\, ■ ■ ■ , B n ]^ = det ([A.Bi, ■ ■ ■ , A.B n ]^j = 

= Yi (-ir + *det([A.B jl ,.-.,ABjj) = 

l<jl<—<jm<n 

E (— l) r+s det(A). det([Sj l5 • • • ,B jm ]) = 

= det (A) E (— l) r+s det^[5j 1 , ■ ■ ■ ,Bj m ]^ = det(A). det(5). 


4. DODGSON ALGORITHM FOR NON-SQUARE MATRICES 

Here, we are in a position to prove the main theorem of this paper by using the previous 
results. To this end, we need to introduce the following notation: 

Notation 1. The (n — k) x (n — l) matrix obtained from A by removing the if 1 , iff, ■ ■ ■ , iff 
rows and the j\ h , jiff , ■ ■ ■ , jf h columns is denoted by (a,;j) ,% k ■ 

■i' ' .Jfc 
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The well-known Dodgson’s algorithm for a square matrix A 
and [15]): 


(a* j) i <i< n , is as follows (see [4] 

1 <j<n 


det 


■ det 

7 ) i^k,k' 


1 <j<n 


3^1,1' _ 
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-h-'H- 

•<ja ^ 
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det 

(^7,7 ) i^k' 

det 

) i^y 
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&i f _ 



(2) 


for all k, k', l, V G {1, 2, • • • , n} with k / k’ and l / V . 

Here, we will generalize Dodgson Condensation Formula for the determinant of non-square 
matrices. 

Theorem 4.1 [Dodgson Condensation Formula for Non-Square Determinants] Let A = (a%,j) 
be an m by n matrix, 2 < m < n, then: 


det 

) 

. det 


= det 

[Cli.j ) i^m 

. det 

(Cli j ) i^m — 1 

l<i<m 


1 

j^n—l,n 

l<i<m 

3=F n 

l<i<m 

j^n-1 

1 <j <n 


1 <j <n 

L J 

1 <j<n 


1 <j<n 



j y£n—l 

. det 

l<i<m 

( a i,j } 

1 i^m — 1 
j^n 

+ det 

l<i<m 

vfim— l,m 

. det 

l<i<m 

(p'iyj) j^n— l,n 


1 <j <n 
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l<j<n 


1 <n 



Proof. We distinguish between two cases: 

Case I. Assume that m + n is odd, and 2 < m < n. Now, we proceed by induction on n. 
Considering Corollary 3.4, for k = 1, we have, 


det 

1 <i<m 
1 <j <n 


(ai,j) = det [Ai, ■ ■ ■ , A n \ = det [Ai - A 2 , ■ 


, A i - A n \ = det 

1 <i<m 
1 <j<n 


(°1 ,j 1 


and, 


det 

1 <i<m 
1 <j<n 



jj^n — l,n 


det 

1 <i<m 
1 <j <n 



jj£l,n — l,n 


therefore, we conclude 


D = det 

( a *b) 

. det 

l<i<m 
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j^n — l,n 


= det 

1 <i<m 
1 <j <n 


( a l ,3 1 


det 
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1 <j <n 


( ^ 1 , jf ^ i^m— l,m 

j 7 ^ 1 , n — 1 , n 


Now, by induction hypothesis 


D = det 

l<i<m 
1 <j <n L 


(Olj fljj) 
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— det 

1 <i<m 
l<j<n 

+ det 

1 <i<m 
l<j<n 


(®ij fljj) 


(d-l j aij'ji^m— l.n 


det 

l<i leqm 
1 <j <n L 

. det 

l<i<m 
1 <j<n 

. det 

l<i<m 

l<j<n 


(ttlj Ojj) i^m-l 


(dlj Oij) i^m — 1 

■77^1, n 


( a bi a *j)jyi,n-l,r 


Here, by simplifying the above expression, and using Theorem 3.3 in reverse direction, the proof 
is completed. 

Case II. We assume that m + n is even. By applying Lemma 3.2, we get 


det 

1 

1 <j <n 


(Pi.j 


det [Ai,-- - , A n \ = det [Ai, • • • , A n , 0 m ] = 


— det [Ai, , A n _ 2 , 0 m , A n _i, A n ] . 


Since the last determinant is determinant of a matrix of size m x (n + 1) and m + n + 1 is odd, 
the result follows from the case I. 
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Corollary 4.2. 

matrix, then: 


[Dodgson Condensation Formula] Suppose A 


det 


. det 



l<i 5< 7 <n 


(oij) be an n by n square 


det 


• dot atopj^n—l] 

— det 

Kp'i.j ) 

. det 

(Olj) i^n - 1 




l<i ,j<n 

j^n-l _ 




Next, we give a generalization of Theorem 4.1. 

Theorem 4.3. Let A = (aij) be an m by n matrix, 2 < m < n,, then: 


det 
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1 <j<n 
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i,3) 
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1<j <n 


{ a i,j 


I ,s f 
j^r,s 


= det 

1 <i<m 
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1 <j <n 


j¥= r 
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( a i,j) if a’ 
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det 

1 <i<m 
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+ det 

l<i<m 

,s' 

det 

l<i<m 


1 <j <n 


1 <j <n 

L J ' J 

1 <j <n 


1 <j <n 



where r,s £ {1, • • • , n},r ' , s' G {1, • • ■ , m}, r j - s,r' ^ s' . 

Proof. In the proof of Theorem 4.1, instead of column one, we choose column k such that 
k r,s. The rest of proof can be treated in the same way. 

Dodgson’s algorithm explains a recursive method for determinant computation. Here, we 
demonstrate an efficient iterative algorithm for computing determinant of non-square matrices 
based on Dodgson’s algorithm. Before describing the proposed algorithm, first we review the 
simplest recursive algorithm which follows directly from Dodgson’s method. The pseudo code 
of the algorithm was shown in Figure 1. 


int Dodgson-Recursive (Mat , m, n) 
inputs : 

Mat: a non-square matrix, 

n: Number of rows in Mat matrix, 

m: number of columns in Mat matrix, 

outputs : 

Determinant of Mat matrix. 

Begin 

if (m == 1) 

Det =0, Sign=l 
for i = 1 to n do 

Det = Det + Sign*Mat (1 , i) , Sign=-Sign 
End for 

else 

DetM(l) = Dodgson-Recursive (Mat (1 :m-l , l:n-l), m-1, n-1) ; 
DetM(2) = Dodgson-Recursive (Mat (2 :m, 2:n), m-1, n-1); 

DetM(3) = Dodgson-Recursive (Mat (1 :m-l , 2:n), m-1, n-1); 

DetM(4) = Dodgson-Recursive (Mat (2 :m-l , 2:n-l), m-1, n-1); 
DetM(5) = Dodgson-Recursive (Mat (1 :m, 2:n-l),m, n-2) ; 

DetM(6) = Dodgson-Recursive (Mat (2 :m-l , l:n),m-2, n) ; 

DetM(7) = Dodgson-Recursive (Mat (2 :m-l , 2:n-l),m-2, n-2); 

Det = DetM(l)*DetM(2)-DetM(3)*DetM(4)+DetM(5)*DetM(6)/DetM(7) ; 

end 

returen(Det) ; 

End 


Figure 1: The simplest version of Non-square Dodgson’s algorithm. 
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Although this algorithm was easy to create and understandable, it is extremely inefficient. For 
time complexity analysis, we will consider the comparison instruction as a basic-operation, and nr 
and n, the size of matrix, as input size. If T(m, n ) be the time complexity of Dodgson- Recursive 
algorithm, then: 


1 


if m = 1, 


T(m, n) 


4 T(m — 1, n — 1) + T{m — 2, n — 2)+ 

T(m, n — 2) + T(m — 2, n) + 1 if n>m> 2. 


It is obvious that this recursive relation has exponential solution. Consequently, the time com- 
plexity of Dodgson- Recursive algorithm is very inefficient. A more efficient algorithm with the 
time cost of polynomial order is developed next, using the dynamic programming technique. 
For l<i,i + l<m,l<j,j + k<n and /, k > 0, let Det(i, i + l,j,j + k) be the determinant of 
matrix Mat = [m Ptq ] i< P <i+i . The pseudo code of the dynamic programming algorithm is shown 

j<q<j+k 

in Figure 2. 


int Dodgson-Dynamic-Programming(Ma, m, n) 
inputs : 

Mat: a square matrix, 

n: Number of rows in Mat matrix, 

m: number of columns in Mat matrix, 

outputs : 

Determinant of Mat matrix 

Begin 

Compute directly all 1 X i determinants with 1 < i < n. 
for 1=1 to m do 
for i=l to m-1 do 
for k=l to n do 
for j=l to j-k do 

Det (i , i+1 , j , j+k) = [Det (i , i+1-1 , j , j+k-1) * Det (i+1 , i+1 , j+1 , j+k) - 
Det (i , i+1-1 , j+1 , j+k) * Det (i+1 , i+j , j , j +1—1) 

+Det(i,i+1, j , j+k-2) *Det (i , i+1-2 , j , j+k)] / Det (i , i+1 , j , j+k) ; 

End 


Figure 2: Dodgson’s algorithm using dynamic programing design. 

According to this pseudo code, it can be shown simply that the time complexity of the proposed 
algorithm is of order of 0((mn) 2 ). Thus, the proposed algorithms has polynomial order time 
cost. This result confirms the importance of the proposed algorithm. 


5. Binet-Cauchy and Trahan formulas for non-square matrices 

In this section, we will extend two well-known determinant identities for non-square matrices, 
Binet-Cauchy and Trahan formulas (see [7] and [14]). 

First, we need the following lemma. 

Lemma 5.1. Let A be an m x n matrix with m < n, then: 

det(A) = ^ a X]x . . . a mjm det 
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where djk is an 1 x n column vector with entry jk equals 1 and the others equals 0. 
Proof. Let A = (a t j), then: 


A = 


E'i=i a iji^ji 

Ylj m =l a mjm^jm 


Since det (A) is linear with respect to each row, 


n 

n 

l 

P 

E 

i 

n 

fiji 

det (A) = E • 

■ ■ E det 


= E a Ll • • • a mjm d e t 


ji=i 

jm = 1 

a mjm*jm 

ilvjm-l 

fijm 


Note that det 


An 


= 0, unless ji, , j m are distinct integer numbers. Furthermore when 


j i, . . . , j m are distinct, obviously det 


An 


= ± 1 . 


L J 

Theorem 5.2. [Binet-Cauchy Formula for Non-Square Matrices] Let A and B be m x p and 
p x n matrices respectively, and m < n, m < p, then: 


det (AB)= J 2 det(T fcl) ... fem )det(s fcl ’-’ fc -), 

i<fci<...<fc m <p 


where det yAki,...kmJ is determinant obtained from the columns of A whose numbers are k \, . . . , k m 

and det is determinant obtained from the rows of B whose numbers are k \, . . . , k m . 
Proof. Suppose that D = AB, that d%j = a^b^j. According to Lemma 5, we have 


det(D) = E dij ... d mjm det 

i < j i < - - . < j m <p 


An 


E (E aikibk 

l<jl<-..<jm<P fcl = l 


Ul 


^ a mkm bkm jm j det 


km — 1 


All 


Aim 


El (°lfel • • • a mk m ) 'y ] (bkljl ■ ■ ■ bkmjm ) ^t 

fel , • . •fkm — i i<ji <...<j m <p 


An 


Aim 


V 

E («ifc 1 ---Omfc m )det(s fcl *-’ fc ™), 


kl — 1 


if ki,.. . ,k m not be distinct, detyB kl, ‘“ ,krn j = 0, then in summation we consider only distinct 
ki , . . . , k m . Since det(B T ( kl b-’ T ( k ™')'\ = sign(r) det(E fcl ’--Am'\ , then for each permutation r of 
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ki , . . . , k m with k\ < ... < k m . we have 

p 

det (D) = E (ai kl . . .a mkrn )det(^B kl ’-’ km ^j = 

kl , ...,/Cm — 1 

= E ( E si S n ( T ) a lr(l) • • ' «mr(m) ) det (. B kl ’-’ km ) = 

l<fcl <...<km<P \ T / 

E det(A fcl) ... )fcm )det(s^--,fc m ). 

l<fcl <...<km<P 


Theorem 5.3. [Trahan Formula for Non-Square Matrices] Assume A = ((Hj) is an m x n 
matrix, where 1 < m < n. Let 


El a i,j > — El > (j — 1 u 


, n 


i= 1 


2=1 


where 


then: 


A rs = (— l) r+s det 

l<i<m 
l<j<n L 




1 n 

det [(ajo)l = — > So-So. 
l<i<m LV m Z-J 3 3 


1 <j<r 


3 = 1 


Proof. The proof is proceed by induction on n. We distinguish between two cases: 

Case I. Assume that m + n is odd. Considering Theorem 3 with X = —A n , and Lemma 3.2, 
we have 

D — det [(Ojj)] — det (Ojj Gin) j^n • 

1 <i<m 1 <i<m |_ 

l<j<n l<j<n 

Now, by induction hypothesis 

n— 1 


^ = -E4E 

m J J 


3= 1 


where 


w 

II 

^i,n) — 

E 

^ ^ ^2,n — 

2=1 


2— 1 

2=1 

ii 

«w 

m 

Y(-l) i+j det 

^ ' l<r<m 

(&r,s &r,n) ■ 

2=1 

2=1 

l<s<n 


m 

= V(-l) i+s det 

^ ' l<r<m 

2=1 l<s<n 

ji 

"c© 

If 

i i 

m 

= E = 

2=1 


s^j,n 


therefore 


1 n— 1 1 n— 1 n— 1 

d = -E(^-^ = -E^--E^ 

m ' m ' m 

i= 1 i =1 i =1 


1 


m 


Y,°i s i- s fX s >■ 


3 = 1 


j=l 
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Now, it is sufficient to prove that Yl'j = 1 = 


•j= 
n m 


= ZB-r* 

j = 1 j — 1 2=1 l<s<n 




EEh) w ,s 

i=i j=i 


1 <r <m 
l<s<n L 


(®r,s) 


= 0. 


Using Lemma 3.1, we get 


l=i 





a i.i 

a 1,2 

* * * &l,n 

det 

(^r,s) r i li 

_ 

= det 

Oi— 1,1 

1 

1,2 

1 

1 

1 <r <m 
l<s<ra 


Cli+ 1,1 

Cti+1,2 

&z+l, n 




. 0"m ,1 

®m, 2 

CLm,n - 


= 0. 


Case II. Now, we assume that m + n is even. By applying Lemma 3.2, we get 

det [Ai, • • • , A n \ = det [Ai, • • ■ , A n , 0 m ] , 
then m + n is odd and the result follows from the Case I. 


6. Conclusion 

An efficient and effective determinant definition for non-square matrices is of great interest 
in linear algebra and applied statistic applications. In this paper, we focus on Radio’s definition 
for the determinant of non-square matrices, because, it has some of the major properties of the 
determinants of square matrices . Unlike previous works, we extract some useful properties of 
this determinant in order to generalize several determinantal identities such as Dodgson, Cauchy- 
Binet and Trahan, for non-square matrices. Also, we present an efficient and straightforward 
algorithm for computing Radic’s determinant based on the generalized Dodgson algorithm and 
a dynamic programing technique that has polynomial complexity. 

More work for defining eigenvalue and eigenvector for non-square matrices and extending 
some matrix decompositions based on these results will be reported in the near future. 
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